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Abstract 

A gauge invariant quantization in a closed integral form is developed over a lin- 
ear phase space endowed with an inhomogeneous Faraday electromagnetic tensor. 
An analog of the Groenewold product formula (corresponding to Weyl ordering) is 
obtained via a membrane magnetic area, and extended to the product of N sym- 
bols. The problem of ordering in quantization is related to different configurations 
of membranes: a choice of configuration determines a phase factor that fixes the 
ordering and controls a symplectic groupoid structure on the secondary phase space. 
A gauge invariant solution of the quantum evolution problem for a charged particle 
in an electromagnetic field is represented in an exact continual form and in the 
semiclassical approximation via the area of dynamical membranes. 

1 Introduction and Overview 

The works by Berezin [|]], Berry ||, [J, and Marinov Q have introduced into mathematical 
physics elegant formulas representing three primary quantum objects (the Weyl non- 
commutative product, semiclassical eigenf unctions, and the evolution wave functions) in 
terms of symplectic area of simple two-dimensional surfaces (membranes) whose boundary 
consists of line segments and pieces of Hamiltonian trajectories in phase space. The area 
of these membranes is determined with respect to the canonical 2-form 

uj = ^J jk dx k A dx j , x = (q,p) G R 2n , J = 

Attempts to generalize some of these formulas to phase spaces with a generic sym- 
plectic form have been undertaken. For symmetric spaces (see M) the opportunity to 
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changes in the geometrical picture. In particular, one can take J to be the matrix 



represent the quantum product via the area of triangle membranes was mentioned by 
Berezin, and it was actually proved in ||, that the semiclassical phase of the product- 
generating kernel is given by such an area in this case. Over Kahlerian manifolds formulas 
for the Wick product and for the solutions of stationary or evolution problems via the area 
of some membranes in the complexification of phase space were obtained in ][7|, || |§ (see 
also flip] ). However, formulas which use only the usual symplectic area are still unavailable 
for the general case. 

In the present paper we analyze three problems related to this topic. 
A first and trivial remark, which one can make regarding generalization of Berezin- 
Berry-Marinov geometrical picture, is that the specific matrix J in the definition of the 
symplectic form u; can be replaced by an arbitrary skew-symmetric matrix without any 

F I 
-I 

where the constant block F represents a homogenous electromagnetic field. But the next 
and more interesting generalization to consider is that of inhomogeneous (not constant) 
tensors F. 

The second natural question is about ordering. All treatments of Berezin, Berry and 
Marinov picture were made for one specific ordering choice: for the Weyl symmetriza- 
tion of the noncommutative coordinates (operators). What happens with other possible 
orderings? 

The third question which our paper addresses is the application of such Weyl and 
non-Weyl symbolic calculus, in the presence of an inhomogeneous field, to solving the 
Cauchy problem and developing a semiclassical representation. 

We begin with the first question and consider the linear phase space R 2n with the 
following inhomogeneous symplectic form 

u F = u + F, F = ^F jk (q) dq k A dq j . (1.1) 

In cases n = 3, n = 4 the form (1.1) describes the structure of the phase space for charged 
particles in an electromagnetic field |Tl|, [E|; the additional summand F is the Faraday 
2-form multiplied by the charge coupling constant e/c. For simplicity of notation we 
include this constant in Fj k ; also note that the order of indices j, k in (1.1) is opposite to 



that used in some textbooks [13]. 



In detail one has the following. In the 3-dimensional case with q = (g 1 , q 2 , q 3 ) then 

F jk = - c e kjl B l (j,k = 1,2,3), 

F = - (V(g) dq 2 A dq 3 + B 2 (q) dq 3 A dq l + B 3 (q) dq 1 A dq 2 ) . (1.2a) 

Here B is the magnetic field, and dF = is equivalent to divB = 0. In the 4-dimensional 
case with q = (q°, q 1 , q 2 , q 3 ), q° = ct one has 

F jk = -e kil B\ (j, k = 1,2, 3), F oj = -E 3 (j = 1, 2, 3), 



F = -[B\t,q)dq 2 A dq 3 + B 2 (t, q) dq 3 A dq 1 + B 3 (t, q) dq 1 A dq 2 ) + eEj (t, q)dq j Adt. ( 1 .2b) 
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The electric field is denoted by E, and dF = is equivalent to the pair of Maxwell 
equations c^dB/dt + curlE = 0, divE = 0; see 0, 0. 

The form up is called the magnetic symplectic form. We show how this form generates 
the Weyl-type associative product * of functions over the phase space. The q 3 ,p k coordi- 
nates on this space correspond to the position of the charged particle and its gauge invari- 
ant kinetic momentum. The commutation relations between the corresponding quantum 
operators q 3 , p k are the following 

[<f , q k ] = 0, [q>,p k ] = ihSi [p v p k ] = ihF kj (q). (1.3) 

The usual realization of these operators in the Hilbert space L 2 (Rq) is q = q, p = 
—ihd/dq — (e/c)$(g), where (e/c)d($dq) = F, and $ is the gauge potential. Specifically, 
for n = 3: 

$ = [Ax, A 2 , A 3 ), curU = B, V j = - e -A 3 {q) (j = 1, 2, 3); (1.4a) 

for n = 4: 

$ = (-a, A h A 2 ,A 3 ), cm\A = B, -{^-— + —j=E, 

Po = -^^ + ^a(t,q), pj = -ih^ - e -Aj(t,q) (j = 1,2,3). (1.4b) 

All the formulae obtained in the paper depend on the symplectic form (1.1) only, but not 
on the choice of potentials, and so, all results are gauge independent. 

The noncommutative product * which we construct necessarily reproduces the com- 
mutation relations (1.3) on coordinate functions: 

q 3 *q k - q k -kq> = 0, q° *p k - q 3 *p k = iML Pj*p k - p k *Pj = ihF kj (q). 

F F F F F F 

The operators of left multiplication q* and p* (as well as right multiplication) are derived 

F F 

without difficulty from these relations, using the noncommutative calculus [|15|, |I6| ; the 
result incorporates Valatin's [17 primitive of the closed 2-form F (see Section 2). We note 



that the problem of a finding a gauge invariant symbol product in a convenient closed 
form was first addressed by Stratonovich |TS|. 



In Section 3 we derive general representations for * in two equivalent forms, both valid 

F 



when F^O. The first is analogous to the exponential Janus derivative representation [IS 



20, |2T[ due to Groenewold. The second is a modification of the (F = 0) Berezin's integral 
for the Weyl product. In both versions there appears an additional electromagnetic action 
(flux) over triangles in phase space. Moreover, the Groenewold-like product formula 
admits generalization for iV multipliers via the magnetic area of polygon membranes. 

Formulas for the * product can be easily expanded to obtain formal h — > power se- 
ries. Higher order terms beyond the Poisson bracket contribution are functions of deriva- 



tives of Fj k . This series coincides structurally with that recently obtained p2| by Miiller. 
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The integral formula for * is also related to constructions |23|, of the Wigner 

F 

function in the presence of electromagnetic fields. In this way in Section 4 we show that 
the -k product can be produced by a convolution over TR n . This convolution is gener- 
ated by a version of the Connes' tangential groupoid |2^, 27, 28 1 but with an additional 
rapidly oscillating factor represented by the electromagnetic flux. Because of the rapid 
oscillations, this type of star-product is outside the framework of the formal deformation 
quantization method. 



In Section 5 we analyze, following the general approach of |[29|| , the structure of the 
symplectic groupoid corresponding to * . This structure is given on the secondary cotan- 
gent bundle T*(T*R n ) = T*R n ®R 2n . The first cotangent bundle, T*R n = is the 
primary phase space, over which we construct the product * . We show how the symplec- 
tic groupoid structure on T*(T*R n ) senses the magnetic correction F in the symplectic 
form, and how the space R 2n , dual to T*IR n , is equipped with a pseudogroup structure 
controlled by the Lorentz momentum of membranes. 

From this point of view, we claim that in the formula for * it would be more natural 

F 

to consider not the usual geodesic triangle, but the triangle with three additional "wings" 
directed vertically (i.e., parallel to the p-direction) in the phase space. The shape of wings 
is determined by the symplectic groupoid structure. 

Then in Section 6 we investigate what happens if the Weyl ordering of noncommuting 
coordinates is changed to some other ordering. The wide (matrix) family of orderings 
introduced in |3(| we relate to phases in the exponential representation of the ^-product. 
These phases can be again presented as symplectic areas of membranes. The membranes 
are combinations of the basic triangle with additional wings that are now not necessarily 
vertical. The shape and direction of the wings exactly control the choice of the ordering 
in quantization, and again it is related to the symplectic groupoid structure over the 
secondary phase space T*(T*R n ). 

The symplectic area of the wings give three additional contributions to the phase. 
The transform from the original Berezin phase (the Weyl case) to the new one, generated 
by the wings, can be considered as a type of gauge transformation of the "symplectic 
potential". On the level of ^-products this is the transformation from the distinguished 
Weyl choice to other ordering choices. In a sense this is the "gauge" of quantization. 

The special features of the Weyl quantization which have made it the preferred choice 
@, [2(], |33], for physical applications are: 1) it treats q and p symmetrically; 2) self- 
adjoint operators have real symbols; and, 3) the Groenewold-Moyal bracket is an even 
function of H, in particular its leading semiclassical correction is 0(h 2 ), not 0(h). From the 
symplectic point of view, the Weyl ordering seems distinguished since the corresponding 
membranes are of the simplest shape (no wings). 

2 1 12 

The Wick normal and anti- normal orderings z* , z and z* , z (where z = q + ip) corre- 
spond to pure imaginary wings of membranes in the product formulas. 

2 1 12 

Other convenient orderings — the standard q,p and anti-standard q,p — correspond 
to the case when the wings are parallel to the basic triangle and the total membrane 
becomes a plane rectangle. These standard and anti-standard cases correspond to the 
push- and pull-groupoid structure on TR n (see Section 4). From the symplectic point 
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of view these cases are singular because of the totally caustic character of the graph of 
symplectic groupoid multiplication corresponding to these cases. 

In the last sections 7 and 8 we apply these ideas to the quantum dynamical problem: a 
charged particle in an electromagnetic field. The basic results for this system were found 
by Dirac, Fock, Peierls, Schwinger. Gauge invariant versions of the WKB approximation 
were developed |31], [34], [5^, |3(| |37|, [38| mostly in terms of integral kernels (Green functions). 



In the context of the present paper we can use the * symbol calculus to obtain a phase 

space gauge invariant treatment of this problem. 

We first study, in Section 7, the pure magnetic situation without electric field. We 
represent the gauge invariant version of the quantum evolution equation over phase space 
and show how the Marinov phase and summation rule are generalized in the presence of 
the magnetic field. Also we investigate how the formulae for semiclassical solutions sense 
the generic wings of the membranes (i.e., arbitrary, not Weyl ordering). 

At the end of Section 7, using the membrane generalization of the Groenewold formula, 
we represent the symbol of the evolution operator exactly in a continual form. This 
continual membrane formula is dual to the Feynman path integral representation. 

Then, in Section 8, we consider a time-dependent electromagnetic field and represent 
the gauge invariant quantum equations over phase space in both the nonrelativistic and 
relativistic cases. Here we use dynamical quantum products which are time-dependent. 
The evolution of commutation relations in time is controlled by the electric field. 

We describe the semiclassical solution of the Cauchy problem using membranes in 7- 
dimensional contact space R ( x x 1^. The boundary of these dynamic membranes are 
given by the solution of two classical systems: one for the given particle and an additional 
one for a "virtual" particle of infinite mass. 



2 Magnetic product for Weyl ordering 

We begin with the definition of magnetic product on a function space over M 2n . The logic 
is the following: we transform relations (1.3) to the standard Heisenberg relations, apply 
the standard Weyl operators of regular representation, and then transform back to the 
magnetic variables. As the result, we obtain a formula for magnetic product in terms of 
left and right regular representation of the algebra (1.3). 
First we recall the properties of closed forms on M. n . 

Lemma 1 Let F = \Fjk(q)dq k A dq> be a closed 2-form on M n . Consider the vector- 
valued 1-form Fdq with components (F(q)dq)j = Fjk(q) dq k and define the two point 
vector potential 

A(q, q') = v -^r f \q - q'\F(q)dq, (2.1) 

\q - q'\ J q > 

where the integral is taken along the straight line path from q' to q, and \ ■ \ denotes the 
Euclidean norm on W 1 . Then for an arbitrary fixed q' G M. n the 1-form A(q,q')dq is a 
primitive of F: 

d(A dq) = F. 
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The choice of a primitive (2.1) (gauge choice) is uniquely characterized by the orthogonality 
condition 



A(q,q')-(q-q')=0. (2.2) 

Note that construction (2.1) is different from that usually used in proofs of the known 
Poincare lemma in the theory of differential forms. On the other hand, (2.1) is just a 
simple particular case of the solution of Lie system related to a general Poisson bracket; 
(see [|16| page 81 and references therein). Formula (2.1) was obtained by Valatin [[[7]] for 
electromagnetic tensors. The characteristic condition (2.2) was stressed by Dirac (see 
p. 101). 

We need several other properties of Valatin's primitive. 
Lemma 2 The following formulas hold: 

d q ,(A(q,q') dq)=d q (A(q',q) dq'), 

A(q,q')-A(q',q)= f F(q)dq, 



A(q, q') + A{q', q) = — ^- f \q - ?\(F(q) - F(q*)) dq, 

\Q <? I Jq m 

where the integrals are taken along the straight line paths, and q* = 2q m — q is the point 
symmetric to q with respect to the middle point q m = |(g + q'). 

Let A be the triangle in M. n with vertices q, q', q" . Consider the integral (flux) of the 
form F: 



F\nx q/I (q } q')= / F. (2.3) 

J A 

Note that here and everywhere in the sequel the orientation of a membrane corresponds 
to the sequence of its vertices (or sides) read from right to left; so the orientation of A 
corresponds to the sequence q <— q' <— q" '. 

Lemma 3 The following formulas hold: 

Flux ,„(?,<?')= [ 9 A(q,q")dq, (2.4) 

Jq> 

d 

— Flux^g, q ) = A(q, q") - A(q, q), 

where the integral in (2.4) is taken along the straight line path. 

Now let us fix the second argument of A(q, q') at some point, say, q' = and introduce 
the operators 

p' =p + A(q, 0). 
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Since q, p satisfy relations (1.3), the new set of operators q, p' satisfy the standard 
Heisenberg commutation relations 



ihd J k , 



W,q k ] = [p' J ,P' k ] = o. 



(2.5) 



Any Weyl-symmetrized function of operators q, p can be transformed to a function of 
operators q, p' by the formula 



where 



f(q,p) = f 



A(q,q') 



3 1 

q ~\~ q ^ ~3i 
— — , p -A(q,q) 



A(qfi + q'(l — /i), 0)dfi. 



(2.6) 



We have used here formulas of noncommutative analysis [[HJ and [|T6j pp. 277-295; the 
superscripts on top of operators denote the order of application. 

Throughout the paper we will not give a characterization of the spaces the symbols 
must belong to in order that (2.6) and the subsequent product formulas are well defined. 
This is a separate technical (and often not simple) question which has been extensively 
investigated in the pseudodifferential operator literature [39 , [10 , 30 , 16| . The reader 
can consider all formulas as formally algebraic or, depending on the formula, assume an 
appropriate simple symbol class such as polynomials, smooth rapidly decreasing functions, 
etc. 

For the Heisenberg algebra (2.5) the operators of left and right regular representation 
are well known. Namely, consider arbitrary Weyl-symmetrized function g' with operator 
arguments q,p': 



g' = g'(q,P') = g' 



q + q t, 
— p 



qg' 




L ' q 


= q + \ihd p i] 


p'g' 






= p' - \ihB q ] 


g'q 


= R^g', 


R' q 


= q - \ihdpi] 


g'p' 




R' p 


= p' + \ihd q . 



(2.7) 



Then the following left and right multiplication formulas hold [g^, [16|, |3(], |41 



(2.8) 



Here we denote d q = d/dq and dp> = 3/ dp'. Operators L',R' satisfying multiplication 
formulas (2.8) are called the left and right regular representation of the given algebra, in 
our case, the algebra (2.5). 

13 

On the right-hand side of the formula (2.6) operators q and q in arguments of A can 
be considered as multiplication by q from the left and from the right, and so they can be 



7 



replaced by L' q and R' q acting on arguments of /. Thus we obtain from (2.6) 

3 1 

/(W) = /'(^P 2 ')=/'(W'), 

where 

f'(q,p') = f(q,p' ~ MKK)) = exp{-A(L' q ,R' q )d pl }f(q,p'), 

and the left arrows mean that operators act on arguments standing to their left. After 
substitution of explicit formulas for L' R' q from (2.8) we conclude 

f'(q,p) = exp{-i(g + \ihd p >, q - \ihd p ,)d p ,}f(q,p). 
Note that by Lemma 2 

pq+u/2 

A(q + u/2,q-u/2)u = / A{q, 0)dq = Flux (g + u/2, q - u/2). 

J q-u/2 

So we obtain the transformation formula 

Proposition 1 Any Weyl- symmetrized function f in operators q,p, satisfying commuta- 
tion relations (1.3) , can be transformed to the Weyl- symmetrized function f in operators 
q,p' satisfying Heisenberg relations (2.5). This transform is given by formula 

f' = U F f, U F = exp{^Flwc (q q-fyhdp)}, (2.9) 
where the Flux was defined in (2.3). 

Using this transform one easily obtains all objects which are needed for the algebra 
(1.3). For instance, the operators of left and right regular representation for (1.3) are the 
following 

L q = Up 1 ■ L' q ■ Up, R q = Up 1 ■ R' q ■ Uf, 

L p = Up 1 -L' p -U F - A(L q , 0), Rp = Up 1 -R' p -U F - A(R q , 0). 
Applying explicit formulas (2.8) for L',R', (2.9) for U F and using Lemma 3 we get 
L q = q + \ihd p , R q — q- Uhd p , 

(2.10) 

L P = p- \ihd q - A(L q , Rg), R p = p+ \ihd q - A(R q , L q ). 

Lemma 4 Operators L q , L p satisfy commutation relations (1.3), and operators R q , R p 
satisfy the conjugate relations {with opposite signs). Operators L commute with R. 

If we know the regular representation of commutation relations, we know the product 
of symbols (see Appendix 2). 
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Theorem 1 The product of Weyl- symmetrized functions in operators satisfying relations 
(1.3) is given by 

f{q,p)-g{q,p) = {f*g){q,p), ( 2 - n ) 

F 

where * is the associative product of functions over~El 2n defined by 

f*g = l p)9 = g(Rq, R P )f- 

F 

We call * the magnetic product of Weyl type. In particular, we have L a = q* , 

F ' ' F 

L p — p*, and R q = *q, R p = *p. These multiplication operators satisfy relations (1.3) 
and their conjugate companions; left and right multiplications commute with each other. 



3 Exponential formula for magnetic product 

The next useful stage in evaluation of the * -product is to bring it into an exponential 

F 



form. First, we apply Proposition 1 and the usual Groenewold formula |L9| , [39f| known for 
Weyl product of symbols over the Heisenberg algebra, namely, 

/'*</ = fexp { - ^J- 1 ^' = f'(q + Tt/2,p) g\q - V/2,p), (3.1) 

where 

"W = ihd p , ~Tt = ihd p , d = (d q ,d p ). 
From (2.9) we calculate 

U F f *U F g = f exp ji Flux (g + (V + lt)/2, q + ( it - V)/2) 

+ l - Flux (q + (lt- *u)/2, q - (V + it)/ 2) - J' 1 J g. 
Applying the inverse transformation Up 1 we obtain 
U F \U F f*U F g) = fexp^-~Flux (q + (^ + lt)/2,q-( < u+lt)/2) 

+ I Flux (g + (IT + lt)/2, q + (lt- *u)/2) 
n 

+ % - Flux (g + {lt- *u)/2, q - (*u + lt)/2) - ^ J' 1 1\ g. 








Figure 1. 
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The sum of three fluxes which appeared at last exponent can be simplified if we look 
at the geometrical picture. Indeed, the three fluxes represent integrals of the form F over 
three triangles with common vertex 0. They are sides of the tetrahedron. By the Stokes 
theorem these three fluxes together are equal to the flux over the bottom triangle; see 
Fig. 1, in which the vectors U\ and u 2 represent the operators if and *u, respectively. 

The relation 

f*g = UF 1 (U F f*U F g) 

F 

and the representation of its right hand side which we have derived above generate the 
following statement. 

Proposition 2 The magnetic product *, corresponding to commutation relations (1.3), 

F 

can be calculated by the formula 

(f*g){q,p) = f(q,p)exp \^-<p(q,ih%,itic%) + y (%'b% - %d^g(q,p). (3.2) 

Here 



0(g,w 2 ,Mi) = / F, (3.3) 

J Aq(u2,«l) 

and A q (u2,u\) is a membrane in W 1 whose boundary is the triangle constructed by the 
middle point q of one side and by the two other sides U2,u±. The magnetic flux (3.3) can 
also be represented in the form 

= [ d^il dv u 2 F(q + (fx - \)u Y + (v - \)u 2 )u 1 . 
Jo Jo 

Note the differential operators in the exponent (3.2) act only on the arguments of 
target functions / and g, but not on the argument q in the flux 0. 

Also note that operators ihd p , which we substitute into the flux in the exponent 
(3.2), are of order 0(h) over the space of non-oscillating functions as h — > 0. So the flux 
actually is of order 0(h 2 ), and the right hand side of (3.2) can be easily expanded as a 
power series in h. Thus for non-oscillating (as h — > 0) functions /, g 



|7|,M=0 1 

^ U,l Ifll— n ) 



|a|,|/9|=0 

where 

e[m + 1] 



es + m 



(s + l)(m + l) (s + l)(s + m + 2) ' 

and e[l] = if / is even, and e[l\ — 1 if I is odd. Here (•, •) is the Euclidean scalar product 
in IR n and the greek letters are multi-indices. 
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Of course, the first two terms of this expansion are 



f*g = fg- ^{f,g} F + 0(h 2 ), 

F 2 

where {, }p is the Poisson bracket on M 2?l corresponding to the symplectic form cup, i-e., 

{f,g} F = d p fd q g - d q fd p g + (d p f,F(q)d p g). (3.4) 

In the case where F is restricted to be the pure magnetic form (1.2a), then the h~ series 
above is equivalent to the one derived in |22[ . 



The formula (3.2) is still not presented in a completely symplectic manner. We have 
there the flux 0, which is the integral of the form F over the triangle in g-space. The 
other part of the exponent in (3.2) can also be related to the area of a triangle, but in the 
(g,p)-space; see Fig. 2 




Figure 2. 

This triangle T, X {V 2) Vj) is constructed by the middle point x of one side, and by the 
two opposite sides V 2 , Vj_. Its projection onto coincides with the triangle /\ q {u 2l u\). 

Theorem 2 The magnetic product can be represented in the form 



{f*9)i x ) = ex P 



h 



S«(V 2 ,Vi) 



wf > f{x 2 )g{x l ) 



(3.5) 



X1=X2=X 



where up is the magnetic symplectic form (1.1), V = (ihd p ,ihd q ) , and x = (q,p) E 



i2n 



In (3.5) one has to evaluate first the symplectic kjp-area of the membrane T lx (V 2 , Pi), 
and then substitute the operators Vy, V 2 for the vectors Vi,V 2 ; the operator Vi is applied 
to the argument x\ and V 2 to the argument x 2 . 

Actually, formula (3.5) may be immediately generalized to the case of several multipli- 
ers. Consider a membrane in M 2n whose boundary is formed by N vectors Vi, V 2 , ■ ■ ■ , V/v 
as sides of a polygon, and by additional point x as a middle point of the (N + l)th side. 
Denote this membrane by T* x (Vn, • • • , Vi); see Fig. 3. 
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Figure 3. 



Corollary 1 The following formula holds: 

C/jv* ...*fi)(x) =exp<^ / uj f \f N (x N ) ■ ■ (3.6) 



F F \ fi 



^x(V N ,...,V!) 



xi=---=xn=x 



An immediate consequence of (3.5) is an integral formula for f*g. Indeed, in (3.5) 

F 

we have a pseudodifferential operator acting on / and g. By the usual formulas via 
Fourier and inverse Fourier transform we easily calculate this action and obtain an integral 
representation for f*g. The triangle H X (V 2 , V\) in this new formula will be described not 

F 

by the sides V 2 ,Vi, but by the middle points x 2 , x x of those sides. Such a triangle in R 2n , 
constructed by three middle points x, x 2 ,x x , we denote by £(x, x 2 ,Xi). 

Proposition 3 The magnetic product of the Weyl type is given by the integral formula 
(f*9){x)= 7 \ 2n L 2n L 2n ex PlT / u F \f(x 2 )g{x x )dx 2 dx 1 . (3.7) 

Actually, the magnetic part F of the form ujp is integrated in (3.5) and (3.7) over the 
projection of S(x, x 2 , xi) onto g-space. 

In the case F — 0, (3.7) becomes the Berezin formula for the Weyl product, and (3.6) 
becomes equivalent to the iV-factor Weyl products obtained in [^C], ^jj . 

4 Tangential groupoid and magnetic cocycle 

Formula (3.7) for magnetic noncommutative product over IR 2 ™ = T*R n can also be derived 
from the Connes type tangential groupoid structure on TR n |26|, [27], 2?|, equipped with 
an additional magnetic phase factor (cocycle). 

The space TM. n consists of pairs (q,u), where 9 6 1" and u G T q M n . Vector u is 
interpreted as a "replacement" at point q. The first and simplest groupoid structure on 
TW 1 is given by formula (q 2 ,u 2 ) o (q x ,u\) = (qi,u 2 + U\) iff q\ = q 2 . This groupoid is 
commutative; it is called the Galileo groupoid. This structure is shown in Fig. 4. The set 
of units of this groupoid consists of all pairs (q, 0); so, the set of units is the configuration 
space R™ considered as the zero-section in TR n . 
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Q2 = qi = Q 
Figure 4. 

Other geometric combinations of the vectors u±, u 2 and the points q±, q 2 may be assem- 
bled to give noncommutative groupoids. Three simple possibilities are shown in Fig. 5. 




q = q\ 

pull-groupoid 
Figure 5. 




push-group oid 




saddle-groupoid 



The multiplication rule, for instance, in the pull-groupoid is (g 2 , u 2 ) o ui) = (q 1: u 2 + Ui) 
iff qi + U\ = q 2 . For the saddle-groupoid one has 

{q 2 ,u 2 ) o (q u ui) = (q,u 2 + Ui) iff q 1 + = q 2 - u 2 /2. 

In this case, q is the middle point of the third side of the triangle: q = qi+u 2 /2 = q 2 —U\/2. 

All these groupoids are specific cases of a general r-groupoid, where < r < 1. 
Pull-, push-, and saddle-cases correspond to r = 1, r = 0, and r = |, respectively. In 
Fig. 6 one can see the generic case corresponding to some r (| < r < 1); in this case 
Q2 — Qi + TU i + (1 — T )u 2 and q = q\ + (1 — t)u 2 . The set of units for all of these groupoids 
is W n q C TW 1 . 




On any (measurable) groupoid there is a convolution of distributions [ 16 , [E? ] 



(0 2 00i)(a) 



(/> 2 (6)0i(c) (iyu 6 (c) 



(4.1) 



a=boc 
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where d/j,b is the Haar measure on fibres of the left groupoid mapping c-*coc 1 . For 
example, in the case of saddle-groupoid we have 

(4>2&(f>i)(q,u) = / (j) 2 (q + u 1 /2,u-ui)(f) 1 (q- {u - u x ) /2,u l )du 1 . (4.2) 



Each convolution over TW n generates a noncommutative product over T*W n just by 
Fourier transform between u and p coordinates: 

r^ u ) = (Av, L elUP/h f^P^ d P=(^ fe* uk f(q,hk)dk. (4.3) 

For example, if one takes the saddle-groupoid convolution © (4.2) then the correspond- 
ing product over T*M n is the usual Groenewold-Moyal product ★ (3.1), i.e., 

(f*g)~ = r®9~. (4-4) 



Remark 1 In this approach the quantum convolution is not an ^.-deformation of a 
commutative one. The deformation parameter h appears only in the Fourier transform 
(4.3); this is what makes it possible to include into the quantization scheme functions 
over T*M n that oscillate rapidly as H — > 0. For non-oscillating functions f(q,p) the 
transform (4.3) admits the classical limit as h — > 0; here one obtains the correspondence 
f(q,p) — > f(q, 0)8(u), h — > between the classical commutative algebra of functions over 
the phase space T*M ra and the commutative subalgebra of distributions over TM. n of the 
type <p(q)5(u) concentrated at the zero "replacement" u — 0. Any replacement by a 
vector uo, represented over TIR n as the distribution 5(u — Uq), after the transform inverse 
to (4.3) generates the phase space observable f(q,p) = exp{—iu p/h}. The corresponding 
operator / (if the magnetic tensor F = 0) is just the Heisenberg translation operator 
fip(q) = ip{q — Uq). This is the simplest example of rapidly oscillating functions / which 
one has to include into algebra of observables over T*M n , but which lie outside the formal 
deformation quantization approach; see also [27|, |4"3 |. 



Now let us return to the relation (4.4) between star-product and convolution. The 
question arises: how to get the magnetic product * (3.7) via the groupoid structure on 

F 

TM n ? Let us slightly generalize (4.1) by introducing an exponential phase factor: 



(02000(a) = / e l ^Mb)Mc)dfx b (c). (4.5) 

J a=boc 

Here $ is a groupoid cocycle, i.e., 

$(d, &o C )-$(do b, c) + $(6, c) - $(d, b) = 0, $(6, c) = -$(c _1 , b' 1 ) 
for any b,c,d G TMJ 1 . The cocycle $ is called the coboundary iff 

$(6, c ) = ^(b) + ip(c) - ip(b o c ), ^(r 1 ) = -tp(b). 
If $ in (4.5) is a coboundary, then convolutions (4.5) and (4.1) are actually equivalent. 
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In the case of saddle-groupoid we can take the following "magnetic" cocycle (cobound- 
ary): 



* = \\ F , ( 4 - 6 ) 

n -lA(q,q 2 ,qi) 



where F is the Faraday 2-form over M n , and A(q, q 2 , q\) is the triangle with middle points 
q,q 2 ,qi- In this case the coboundary function ip at the point b = (q,u) is given by the 
integral ip(b) = I J Adq along the chord [q — hu,q + \v\. Note that ip controls the 
additional phase factor in the gauge invariant version of the Wigner function p4|, p5 . 



Proposition 4 Let © be the saddle-groupoid convolution, equipped with the magnetic 

cocycle (4.6). Then the Fourier transform (4.3) relates this convolution to the magnetic 
product -k via 

(f*g)~ = r®g~- 

F F 

Note that other types of groupoids (like pull-, push- ) will also generate some magnetic 
products over M 2n = T*M n . In Section 6 we will consider a variety of such products. On 
the other hand, we will show that not all of products over T*M n are generated from the 
convolution over TR n in this fashion. 

In Section 5 we analyze the magnetic product (3.7) from the view point of symplectic 
groupoid structure of the secondary phase space (on the secondary cotangent) T*(T*R n ). 
Regarding this we mention the following. 

Proposition 5 Each groupoid structure on TM. n whose set of units is uniquely de- 
termines a symplectic groupoid structure on T*(T*M n ). Here T*M. n is equipped with the 
symplectic form uq = dp A dq. If, in addition, on TMJ 1 a cocycle of type (4.6) is given, 
then on T*(T*W n ) we have a symplectic groupoid structure corresponding to the magnetic 
form ujp. 



5 Symplectic groupoid and membranes with wings 

Formula (3.7) represents the exact magnetic product. On the other hand, it is known 
in a very general context, how to construct quantum products of functions over an 
arbitrary Poisson manifold M in the semiclassical approximation, to all orders in h — > 0. 
This approximate product takes the form 



Cf*#)(z)-/ / K h (x,x 2 ,x 1 )f(x 2 )g{x 1 ) dx 1 dx 2 . (5.1) 




Here x,X2,x\ G A/", and K% is a "wave function" corresponding to some Lagrangian 
submanifold A* in the "phase space" £ x £ x S, where £ is the symplectic groupoid over 
J\f. In our case N = R 2n = T*R n with Poisson bracket (3.4). The Lagrangian submanifold 
A* is the graph of groupoid multiplication in £. If this graph is one-to-one projected onto 
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the "configuration" space M x jV x jV along the polarization then the function Kn is just 
a WKB function 

K h = exp{iS/h}(p + 0(h) (5.2) 

whose phase S is the Poincare-Cartan action on A* and (p is the solution of the corre- 
sponding transport equation. 

In specific cases, for instance, in our case Af = R 2n with bracket (3.4), the asymptotic 
formula (5.1) becomes exact and the remainder 0(h) in (5.2) is absent. Indeed, let us 
compare (5.1), (5.2) with formula (3.7) for the magnetic product * . First, one can describe 

F 

the symplectic groupoid £ in our specific case. Let us set £ = T*R 2 ™ = R 2 ™ © R 2n , and 
denote by (x, y) points in £. The space R 2?1 is imbedded into £ as a zero section {y = 0}. 
The symplectic form on £ is dy A dx. The operators L, R (2.10) of left and right regular 
representation of algebra (1.3) generate two mappings 

I : £ -> R 2n , r:£^ R 2n . (5.3) 

Here I = (l q , l p ) and r = (r q , r p ) are just symbols of L = (L q , L p ) and R = (R q , R p ), i.e., 

L = l(x, —ihd x ), R = r(x, —ihd x ). 

In our case formula (2.10) reads 

l q (x, y)=x q - y p /2, l p (x, y) = x p + y q /2 - A(l q , r q ), 

(5.4) 

r q (x, y) = x q + y p /2, r p (x, y) = x p - y q /2 - A(r q , l q ). 

Mappings (5.3) are Poisson and anti-Poisson, i.e., / preserves brackets, r changes the 
sign of brackets (recall that on R 2 " we have the bracket (3.4), and the bracket on £ 
corresponds to the symplectic form dy A dx). 

The groupoid structure on £ is defined as follows: points m 2 ,mi G £ are called 
multiplicable iff r(m2) = ^(^i); the product m = m 2 o m 1; by definition, is a point in £ 
such that l(m) = l(m 2 ), r(m) = r(mi). The subspace R 2 ™ C £ is the set of units of this 
groupoid, and mappings (5.3) are left and right reduction mappings: 

l(m) = mo m , r(m) = m~ l o m. 

The graph A* C £ x £ x £ of this groupoid multiplication consists of all multiplicable 
points and their products 

A* = {(m, 777.2, mi)\m = m 2 o mi}. (5.5) 

If on £ x £ x £ we introduce the symplectic form dy A dx — dy 2 A dx 2 — dy\ A dx±, then 
the sub manifold A* is Lagrangian (see details in [^, |44}j). 

In the case F = (i.e., A=0) formulas (5.4) are interpreted as "middle point of chord" 
relations: x = \(l + r), y = JV, where V = I — r G TR 2n . So, the groupoid structure on 
R 4n = T*R 2n is given by the triangle rule (see Fig. 7): 

m = m 2 om 1 , m = (x,y), m 1 = (x 1 ,y 1 ), m 2 = (x 2 ,y 2 ); 

(5.6) 

y = yi + V2, x = Xi + \J 1 y 2 = x 2 -\J 1 y 1 . 
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Figure 7. 



For arbitrary given triple of points x,X2,X\ G K. n we uniquely construct the triangle 
for which these points are the middle points of its sides, and so reconstruct elements 
m,m2,mi G £ such that m = m 2 o m\. This means that the graph A* is one-to-one 
projected onto R 2n x R 2n x R 2n along the "vertical" y-polarization. Hence, the kernel K h 
has the WKB form (5.2), and its phase is 



-'(0,0,0) 

Here y,y2,Ui are determined via x, X2,x± following the triangle multiplication rule; the 
initial point (0, 0, 0) corresponds to triple of elements m Q = m Q om , where m = (0, 0) G £. 



The integral (5.7) is taken over an arbitrary path on A* connecting the triple (0, 0, 0) 
(i.e., the degenerate triangle) with the triple (x, X2,xi) (i.e. the given triangle). This path 
is actually a family of triangles in M 2n ; see Fig. 8. The phase (5.7) in the case F = is 
just equal to the area of the final triangle 



Now let us see what happens in the general magnetic case F^0. From (5.4) it follows 
that the groupoid structure now is given by formulas 




(5.7) 




Figure 8. 




\{l + r) + V F , y = JV + Y F , where V = I - r. 



(5.8) 
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Here we have introduced notations 



V F = (0; A 8 ) G M™ © Rp, Y F = (A a ; 0) G E" © R" ' 

(5.9) 

A s = §(A(Z 3 , r 9 ) + A(r„ /,)), A a = A(Z„ r g ) - A(r„ Z ? ). 

Comparing with the multiplication rule (5.6) in the case F = we see that points 
x,X2,xi are no longer the middle points of sides of the basic triangle and do not even 
belong to those sides. They are shifted in the p-direction by the vectors Vp. 

Another difference from the case F = is that the usual multiplication rule y = y\ + 
(see (5.6)) fails to hold for inhomogeneous magnetic case. New "magnetic" rule is 

y = yi + y2 + ([ VF ; o), (5.10) 

J A 

where the vector- valued closed 2- form VF is defined by VF = ^VFj k dq k A dqi , and A = 
A(q,q 2 ,qi) is the triangle in g-space with middle points (q, (72,91), which are projections 
of x, x 2 , xi G R 2n onto Rg. Note that 




Lorentz momentum. 



The vector 2-form VF is a measure of magnetic inhomogeneity. So, this form controls the 
modification of the usual (F = 0) multiplication rule. The space Ry n (dual to E^ n ) now is 
not even a group, it is a pseudogroup over the Poisson manifold R^ n (see details in general 
case in |^6| , the word "pseudo"- reflects the fact that the product (5.10) in the y-space 
depends on ^-coordinates as additional parameters, and the condition of "associativity" 
of the product (5.10) senses this dependence). The non-trivial part of the pseudogroup 
structure (5.10) is determined by the momentum of the membrane A in the magnetic 
field. 

Formula (5.7) in the case F 7^ still represents the phase in the product (3.7) if we 
put there the corrected values y = JV + Yp. But now the triangle no longer represents 
the groupoid multiplication rule, and from this point of view, it seems unnatural to keep 
this triangle in the formula for the phase. 



X2 Xl 




Figure 9. 



Actually, we have another configuration (see Fig. 9) related to the triple (x,X2,xi). 
This configuration consists of triangle S(£, X2,xi), where x = x — Vf, and also of three 
additional triangles directed "vertically" (parallel to p) with "top" vertices x,x 2 ,Xi. We 
call them magnetic wings. 
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A magnetic wing is characterized by a sequence of its vertices [l,x, r] related to each 
other by (5.8) (see Fig. 10, left picture). The configuration of membrane with wings, as 
we see, exactly corresponds to the new groupoid multiplication rule for F^O. 





x 

Figure 10. 



Lemma 5 



/ uop = I up, / up — 0. 

JH{x,X2,xi) JH(x,X2,xi) J vertical wing 

Both statements of the Lemma follow from the orthogonality condition (2.2). 

Denote by Yip(x, x 2 , x±) the triangle S(x, x 2 ,x±) together with three magnetic wings 
described above. The boundary of this figure consists of six straight line segments. The 
figure itself looks like inflected hexagon. More generally, T,p could be any membrane in 
the phase space 1R 2 ™ with that six segment boundary. In view of Lemma 5 we have 



UJp = / iOp. 

H(x,X2,Xl) J T,p{x,X2,X\) 



Theorem 3 The Weyl type magnetic product over T*M™ is given by 

(f*d)(x) = , \ 2n I I 6X p|t / uj F \f(x 2 )g(x 1 )dx 2 dx 1 , (5.11) 

where T*p is the wing membrane corresponding to the magnetic groupoid structure (5.8) 
on T*(T*R n ). 



6 Ordering in quantization 

Now we demonstrate that membrane wings introduced in previous section are actually 
very natural objects in the quantization framework. We show how different configurations 
of wings relate to different choices in the ordering problem. In particular, we'll see why 
the Weyl ordering choice (Weyl symmetrization) looks like an optimal choice. 

Let us take a constant real In x 2n matrix M (actually, a linear operator in the space 
tangent to T*IR n ) which obeys the condition 

M T J+JM = 0, (6.1) 
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where M T is the transposed matrix. One can represent M by its n x n blocks as follows 

f N K 

Now let us take an arbitrary function / = f(q,p) (say, a polynomial) and determine 
the following general ordering of operators q,p' (generators of the Heisenberg algebra): 

JM = f (l N fi 1 ( | + l R{ p _ pr^ J, _i s( i _ | _ i^T^, _ ^ ( g_ 2) 

We would like to obtain the product formula 

f M -g M = k M , k = f M *g. (6.3) 



To derive the product * we first note that the M-ordering (6.2) is related to the Weyl 
ordering via 

f M =(U M fMp'), (6.4) 

where 

U M = expj^S M (y)j, V = (ihd p ,ihd q ), 
and the function S M is defined by the matrix M as follows 



S M (V) = l-(JMV,V)= [ u 
* Ja m (v) 



X\=X2=X 



Here A M (V) is the triangle, called the M-wing, generated by vector V and by vector MV 
applied at the middle point of V (see Fig. 10). 

Using U M we can calculate the product (6.3) by the formula 

f%g = U- M (U»f*U M g), 

where * is the Weyl product (without magnetic correction, at first). Looking at exponen- 
tial representation of the Weyl product, one concludes that 

f*g = e W {^ M (V 2 ,V 1 )}f(x 2 )g(x 1 ) 
where V = (ihd p ,ihd q ). In this formula we use the notation 

® M (V 2 , Vi) = / co + S M (V 1 ) + S M (V 2 ) - S M (V 1 + V 2 ), 

JX(V 2 ,Vi) 

where ^(V 2 , Vi) is the triangle generated by vectors V±, V 2 . 

The phase function $ M can be written as the symplectic area of the wing membrane 
S^(V 2 , Vi) generated by four triangles: 

(V 2 , V,) = A(V 2 , V,) |J A M (V 1 ) |J A M (V 2 ) |J A M (V l + V 2 ). 
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Its boundary consists of six line segments. For reasons of uniformity in notation the 
point x is shown; actually, the uj -aiea of T,^, 1 (V 2 , V\) is independent of x. 

The picture for is the same as in Fig. 9, but wings now are M- wings as in Fig. 10 
(right picture). 

Theorem 4 For an arbitrary 2n x 2n matrix M , satisfying (6.1), there is a star-product 
over M 2n given by the wing membranes: 

if * 9){x) = exp \\ \ uo\f(x 2 )g(xi) . (6.5) 



h 



Xl=X2=X 



This product corresponds to the M -ordering rule (6.2) of non-commutative operators q,p' 
(generators of the Heisenberg algebra), so that formula (6.3) holds. 

The family of orderings (6.2) was introduced and studied in detail in |!0|, where the 
following pseudodifferential formulas for M-product were obtained: 

(/ * 9)(x) =f[x -ih{\ - M)J- 1 i )g{x) =g(x +ih{\ + M) J~ l d x )f(x). 
In the particular case 

M=a-r)(l ° r ) (6.6) 



^ ' \ -I / 
the ordering (6.2) simplifies to 

1 3 2 

f T = f(rq+(l-r)q, p'j. (6.7) 

The family of r-wings corresponding to this family of orderings is represented in Fig. 11. 
The whole membrane for the specific case r = is shown in Fig. 12. 





Figure 11. 



Figure 12. 



The cases r = and r = 1 are called the standard and anti-standard ordering choices 
(they correspond to push- and pull- groupoid structures on TM. n , see Section 4), and the 
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case r = i is the Weyl ordering (corresponding to the saddle-groupoid structure on TM. n ) . 
In the later case M = and wings in membranes S M are absent. 

Matrix M is assumed to be real, but if we formally take M = iJ/2 (i.e., N = 0, 
S = —i/2, K = i/2), then the transform U M = exp{— h{d q +dp) /A} relates Weyl ordering 
to the Wick normal ordering choice. In this case wings are pure imaginary (see Fig. 13), 
and the membrane S M coincides with hexagon membrane introduced in [7|, ||. 




Remark 2 In the general case, block iV of the matrix M controls the position of the 
point g6R n with respect to the middle point of the vector u in the TM n -groupoid inter- 
pretation. The block S controls an additional TM n -groupoid cocycle (in ^-coordinates) 
which appears in the convolution formula (see Section 4). In contrast to that, the block K 
in the matrix M generates the transformation of distributions over TM™ of the following 

M 

type: exp{— ihKd q ■ d q /2}. This is not a point transformation. Thus the product * cor- 
responding to matrix M with K ^ can not be obtained from the groupoid convolution 
over TW" 1 by construction of Section 4. In particular, the Wick product is of such type. 

Now let us consider the symplectic groupoid structure on the secondary phase space 

M 

T*(T*IR n ) corresponding to the * product (6.3). In the same way as in Section 2 we 
calculate for the Heisenberg algebra (2.5) operators of left and right representations cor- 
responding to the ordering choice (6.2). We know the transformation from (6.2) to the 

Weyl ordering; it is given by operator (6.4). Thus the left operators L^ 1 = q * , L^f = p' * 
and right operators — * q, M*/ = * p' are given by 

L M = jj-M . L , . jjM^ R M = V -M . R > . (6 g) 

where L',R' are determined in (2.8). We represent L M ,R M via their symbols l,r: 

L M = l( x , -ihd x ), R M = r(x, -ihd x ), 
and easily calculate /, r by (6.8) and by the definition of U M . The result is 

l(x, y)= x + {\- M).r l y, r(x, y)= x -{\ + MjJ^y. 
The inverse mapping (l,r) (x,y) is given by 

x = \{l + r) + MV, V = l-r = J~ x y. (6.9) 
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For example, for the ordering cases (6.7), formulas (6.9) are represented in Fig. 11. 

The corresponding symplectic groupoid structure on T*(T*R n ) = IR^ n © Ry n is given 
by the rule represented in Fig. 14. We observe that this structure is exactly given by the 
membrane with wings S M which we described at the beginning of this section. What is 
new now is that we have identified the positions of points x,x 2 ,x± exactly as vertices of 
wings of the membrane S M . 




Figure 14. 



Corollary 2 Let det(| — M) ^ 0, then the integral version of the product formula (6.5) 
reads 

{! * 9)W = p^Sj det(i - 501 / / 6XP ^ !<<„,,„ , *}><*>*<*>' fa * fa " 

(6.10) 

where T, M (x, x 2 , xi) is the membrane with three M-wings having vertices x,x 2 ,xi. 

Note that for the r-ordering case (6.6) the denominator in formula (6.10) det(| — M) = 
r n (l — r) n is zero if r = or r = 1. We see that the cases r = (standard ordering) and 
r = 1 ( ant i- standard ordering) are special. In these cases the graph A* of the groupoid 
multiplication (see Section 5) is not one-to-one projected onto the "configuration space" 
M^ n x M. 2 x n 2 x R^™, i.e., the graph is totally caustic and the integral kernel K h (x, x 2 , x±) (5.1) 
is not of WKB-type (5.2). This is one of essential diffferences between representations 
(6.5) and (6.10) for the *-product. 

Corollary 3 For several multipliers there is the following product formula 

(f N *...*f 1 )(x) = exp\l uj \f N (x N ) . . , (6.11) 

k 11 JliM(y N ,...,Vl) > Xi=-=Xjv=X 

where H^(Vn, . . . , Vi) is the membrane with N + 1 wings. 

These formulas are naturally generalized for products of differently quantized multi- 
pliers, i.e., for the case when each multiplier has its own ordering choice: 

(if i ( 6 - 12 ) 

^)=exp{-/_ . ) / ,(^)... /,(,,) 



Xl = ---=X]\[=X 
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Here the membrane y^ n+u '"' Mi is constructed by N + 1 wings each of different configu- 
ration determined by different matrices M N+ i, . . . , M x . 

In conclusion of this section we consider the magnetic case F^O. For simplicity, we 
concentrate on the r-ordering choice 

/ r = /(4+(l-rH I), (6.13) 

where q, p satisfy the magnetic commutation relations (1.3). 

We first transform the magnetic r-ordering to the magnetic Weyl ordering: 

r = (U T f)(q,p), CT = exp{(i-r)«a,}, (6.14) 

where u = ihd p . Then for product of two r-ordered observables we have from (2.1 1): 

t ■ ¥ = (U T f)(q,p) ■ (U T g)(q,p) = (U T f t U T g)(q,p), 

F 

where * is the magnetic product (3.2), (3.5), corresponding to the Weyl ordering. Trans- 

F 

forming back the Weyl symbol to r-symbol by the transform (U T )~ 1 , we obtain the product 
formula 

t ■ 9 T = k T , 



where 



k = f ★ g = (U T )-\U T f*U T g). 

F F 



By formula (3.2) we calculate the new r-magnetic product as follows 

(f*9)(q,p)= (6-15) 

F 

exp{(r- i)ua,}^exp{i0(g,U2,wi) + \{uid q2 - u 2 d qi )} 
x ex P{(| - T )(uid qi +U2d q2 )}f(q 2 ,P2)g(qi,Pi) S ) 
= exp + (r - i)(wi + u 2 ),W2,«i) + \{uid q2 - u 2 d qi ) 

(r- ^)^(u 1 + u 2 )(d qi +d q2 ) -uid qi - U2d q ^)}f(q 2 ,p 2 )g(qi,pi] 



qi=q2=q,pi=P2=p 



+ 



qi=q2=q,pi=P2=p 



Note that the flux <p(q,u 2 ,ui) was defined by (3.3) via the triangle A q (u 2 ,ui) in R" 
with the middle point q of one of sides. Now we see that the position of the middle point 
q differs from q by an additional vector (r — 1/2) (ui + u 2 ). So, 

q = q+{\-T){u l + u 2 )=q + (MV) q , 

where M is matrix (6.6) corresponding to r-ordering, and V — Vi + V 2 is a vector in R 2n 
whose g-component is u = u x + u 2 . 
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The total phase in (6.14) is equal to the magnetic area of the membrane in IR 2 " con- 
structed by the triangle with sides V 1: V 2 , V\ + V 2 and by three additional wings generated 
by matrix M of special type (6.6). The configuration of the wings in this special case is 
shown in Fig. 11. 

The position of the vertex of the wing is x T = x+MV, where x = (Z+r)/2. But actually 
this position should be changed to match the symplectic groupoid multiplication rule. To 
find this rule we have to calculate the operators of left and right regular representation 
in our r-case. In view of (6.14) these operators are given by formulas 

L T = (CT)- 1 • L ■ U T , R T = (U^- 1 ■ R ■ U T , 

where L,R are operators (2.10), corresponding to the magnetic Weyl ordering. From 
(6.14) and (2.10) we obtain 

L q = q + Tihd p , R T q — q — (1 — r)ihd p , 

L; = P -(l-r)ihd q -A(L' r q ,K r q ), R; = P + rihd q -A(K r q ,L r q ). 

So, if we represent these operators by symbols U = l(x, —ihd x ), R T = r(x, —ihd x ), 
then two groupoid mappings appear 

I : T*R 2n -> R 2n , r : T*IR 2n -> R 2n , 

where 

l{x,y)=x + {\- M)J- l y - (0; A(l q , r q )) , 
r(x, y) = x - {\ + M)J-'y - (0; A(r q , l q )) . 

From (6.16) we reconstruct x and y via l,r 

x=±(l + r) + V F + MV, y = JV + Y F , (6.17) 

where V = I — r, and vectors Vf, Yf are given by the same formulas as in (5.9) (i.e., the 
same as in the case r = 1/2 or M = 0). 

Relations (6.17) determine the final configuration of the magnetic r-wing. Position of 
the vertex x of this magnetic wing is shifted by vector Vf with respect to position x T . 

Lemma 6 The magnetic area of the magnetic r-wing with the vertex x is equal to the 
magnetic area of the wing with vertex x T . 

In view of this lemma the phase (6.15) can be represented by the area of a membrane 
S^(V2,Vi) constructed by the triangle of vectors Vi,V2,Vi + V2, and by three magnetic 
r-wings over each of these vectors. The point x is the vertex of the magnetic r-wing over 
Vi + V 2 . 

Theorem 5 The magnetic star-product corresponding to r-ordering (6.13) of noncom- 
mutative coordinates q, p is given by formula 



(/ * 9){x) = exp \ l - / u JF }f{x 2 )g{x 1 ) 



'Sj(V 2 ,Vi) 

where YZ £ is a membrane with magnetic r-wings. 



X\=X2=X 
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The immediate corollaries from this statement are formulas for several multipliers and 
also the integral formula for * via the magnetic area of membranes S r (a;, x 2 , x\) supplied 

F 

with magnetic r-wings. These corollaries are formulated by the same way as (6.10)-(6.12), 
but with magnetic form ujp in the exponent. 

7 Dynamics via membrane area 

Quantum dynamics of a charged (spinless) particle in electromagnetic field can be de- 
scribed 

- in the nonrelativistic case by the Schrodinger equation 



dtp 



1 ( - t d e .\2 

n I A l +ea 

2m \ oq c / 



ij, (7.1) 



in the relativistic case by the Klein-Gordon equation 



(^ft— -eaj i) = (? y — ih— - -A J + mV tj), (7.2) 

where A and a are magnetic and electric potentials of the field, and metric is assumed to 
be Euclidean. 

First, we consider the pure magnetic time- independent situation a — 0, A — A{q). 
Let us introduce operators p = —ihd/dq — -A(q). Then the dynamical equations can be 
reduced to studying the evolution operator exp{—itH(p)/h}, where H{p) = p 2 /2m in the 
nonrelativistic case, and H(p) = ±c^p 2 + m 2 c 2 in the relativistic case. 

More generally, in a presence of an additional non-Euclidean metric (gravitational 
field) the Hamiltonian will depend on the g-coordinate as well: H = H(q,p). Say, H ~ 
gi k (q)pjPk/2m in the nonrelativistic case. 

So, the general problem is to study the operator 

U t = exp{- l ±H}, H = H(q,p). (7.3) 

In particular, we are interested in its asymptotic behavior as H — > 0. 

Note that the symbol H and the permutation relations between quantum coordinates 
q,p (1.3) are independent of the gauge choice of the magnetic potential A. Thus the 
semiclassical approximation for [/* as h — > 0, written in terms of the phase space sym- 
bol H(q,p) and the phase space noncommutative structure (1.3), is automatically gauge 
invariant. 

Let us represent the operator U l in the Weyl form 

U^tf^p), 

then for symbol U l one obtains the following equations 

th^- = H*U\ U° = l. (7.4) 
at f 



d e ,\2 
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Using Theorem 1 we transform (7.4) to a pseudodifferential form: 

ih ^!M. = Hh{Xj -ind x )u\x), u° = 1. (7.5) 

Here TCn(x,y) is the Weyl symbol of the operator H(L q ,L p ), and L q , L p are given by 
(2.10). Obviously 

n h = n + o(h 2 ), n = H(i), (7.6) 

where / is the Weyl symbol of L, i.e., L = l(x, —ihd x ). The explicit formulas for / = (l q , l p ) 
are found in (5.4). 

In view of (7.6) the principal term of the semiclassical solution of the Cauchy problem 
(7.5) is determined by the Hamilton function 7i - For small enough time interval t e [0, T] 
the approximate solution has the simplest WKB-form 

U\x) = exp^S(t,x)^u t (x) + 0(h), (7.7) 
where the phase S is the solution of the Hamilton-Jacobi equation 

§ + «„(*,!) =0, Sl r 0, (7.8) 
and the non-oscillatory amplitude u l is the solution of the "transport" equation 

M + l(^U, B 4V)=0, ^ = 1. (7.9) 



dt dx \ dy V ' dx / 
In order to solve (7.8), (7.9) one has to consider the Hamiltonian system 



<)Y ' ' " dX ' 



= x°, Y 

t=o 



= 0, (7.10) 

t=o 



and determine 

(YX-H )dt, «'=(^o) , (7-11) 

where x° = (q°,p°) is taken from the equation 

X(x°,t) =x. 

The time interval t e [0, T] over which the Jacobian VX/Vx° >5>0 (for any x° e M 2n ) 
is exactly the interval where the solution U l can be represented in the WKB form (7.7). 

Lemma 7 Trajectories of (7.10) are given by formulas 

X(x°,t) = l( 1 t (x )+x )+v t (x°), 

Y(x°,t) = J(7*(x°) -x )+y*(x°), U j 
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where 



0; \ {A{% q°) + A(q°, 7*)) ) , y* = (a( 7 <, g°) - 7*); o) , 



and 7* = (7* 7*) zs £/ie trajectory of the Hamiltonian system corresponding to the function 
H and the magnetic Poisson bracket (3.4) 

dH, s . OH . . w .dH. . _ in . 
lq = ^T (7<? ' 7p) ' 7p = ~~dq~^ q,lp ' ~ F (^)^(7 g ,7 P ), (7.13) 

7 =a; o = 
t=0 

The proof of this lemma follows from (5.8) and from the fact that TCq = H(l), and so 
components of the mapping r : M^ n © R 2n IR 2ri (5.4) are integrals of motion for system 
(7.10), i.e., r(X,Y) = r(x°,0) = x° is constant in time. 

After substitution of (7.12) into (7.11) one obtains the following result. 

Theorem 6 Let q,p satisfy commutation relations (1.3) with magnetic tensor Fkj, and 
H = H(q,p). Then for small enough time t the semiclassical approximation for the 
magnetic Weyl symbol of the evolution operator exp{ — j?H} = W*(g,j3) is given by the 
formula 

J- 1 / 2 exp {LJ u F - + 0(h). (7.14) 

Here the membrane £ = S* (see Fig. 15, left picture) is constructed from the piece of the 
Hamilton trajectory (7.13), which connects points x° and 7 *(a; ), and from the magnetic 
wing with vertices [ 7 *(x°), x, x°], where 

x= I( 7 *(x°)+x°) +v t (x°). (7.15) 

Here 

v\x Q ) = (0 ;^(A(Y(x ) q ,Q ) + A(g°,7V) ff ))). x° = (q°,p°), 

and A is the Valatin primitive (2.1). The Jacobian J in (7.14) is given by 

J(x°,t) = det[|(/ + rf 7 '(a; )) +dv t (x )}, J > 5 > for t e [0,T], 

the Hamilton function H in (7.14) is evaluated on the trajectory, i.e., H = H(x°), and 
x° is assumed to be expressed in terms of x,t via the equation (7.15). 
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Remark 3 The group property of the family of symbols U l over R 2n reads 

U t2 * U tx =U t2+tl . (7.16) 

F 

In terms of WKB-phase functions (7.14), the identity (7.16) requires that 

ujf + / tup + up = I ujf (7-17) 

(see Fig. 16), where Ylp(x, x 2 , X\) is the hexagon membrane with magnetic wings, defined 
at the end of Section 5. Note that in the case F = formula (7.17) coincides with 
the phase addition rule obtained by Marinov Q. In that particular case the magnetic 
"anomaly" v* in (7.15) is absent and the magnetic wings of the membranes disappear. 




Figure 16. 



Also note that the Hamilton function H could be time-dependent. In this case the 
first membrane phase factor in formula (7.14) is the same, but the second phase factor 
becomes exp{ — | J Q * H{pf{x°),t) dt}; the trajectory 7* is now the solution of system (7.13) 
with time-dependent Hamiltonian H. 

Remark 4 One can use not only Weyl but any other ordering choice to represent the 
evolution operator as a function in coordinates q,p. Then formula (7.14) still holds with 
membrane £ constructed by wings corresponding to the given ordering choice (see Sec- 
tion 6); equation (7.15) and the Jacobian J are changed following (6.17). Moreover, in 
||30|1 it was proved that using and combining different orderings it is possible to avoid the 
difficulty with time limitations t 6 [0, T] where the WKB-approximation works. 

For the Wick ordering choice the wings are pure imaginary (see Fig. 13) and the 
membrane representation (7.14) coincides with those obtained in In this case the 
Jacobian J is never zero and representation (7.14) is global in t. 

In conclusion of this section we apply formula (3.6) to derive the symbol U l not 
asymptotically but in an exact continual form. Namely, one can use the approximation 
yt/N _ eX p 1 _ jLh] + 0(N~ 2 ) and obtain the Trotter type formula 

U l = lim exp { — ^—h\* exp ( — t^hX 
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(N multipliers). Applying (3.6), one derives 



ifix) = lim exp { 



X x (V N ,...,Vi_) 



u F > exp < - 



it 

hN 



N 



(7.18) 



XI- 



= Xff=X 



Now the question is how to represent this formula in a continual form. 



Note that each vector field v on M 2n and any point x E 



o2n 



t e 1 determines a 

membrane S* (y) C BL 2n whose boundary is constructed from a piece {T^ | < fi < t} of 
the trajectory of the field v in IR n and from the magnetic wing with vertices [r',x, T ] 
(or, the magnetic r-wing if one wants to use the general r-ordering). The integral over 
the membrane in (7.18) is an approximation of the integral over S^(t>) with a convenient 
choice of v. 

Theorem 7 The following continual formula for the symbol of the evolution operator 
(7.3) holds: 



U t (x)=exp<— / ujp\ exp j — — / H (x(fi)) dfA 



(7.19) 



x(fi)=x 



Here (x(/i) = (q(fi) , p(/^)) | < /i < t} are continuous paths in M 2n andv= («M/5p(/x), 
ih~5/5q(ii)) is the variational derivative operator acting on the path functional. 



Formula (7.19) is dual to the Feynman path-integral formula p5j f46| 



HI. The 



difference between (7.19) and the path integral is the same as between (3.5) and (3.7). 
The known Wick and Hori formulas |49|, [5(| for the symbol of the evolution operator 
(see also generalizations in plf ) are structurally close to (7.19), but use a different first 
exponential factor. The membrane exponential factor in (7.19) clearly demonstrates the 
influence of the magnetic form up to the quantum dynamics. 



8 Electromagnetic fields and space-time membranes 

Let us now consider general time-dependent case, i.e., A = A(t,q), a = a(t,q) in (7.1), 
(7.2). We again study the Cauchy problem for the Schrodinger or Klein-Gordon equations. 

As a first step one can remove the electric potential a from equations by introducing 
the new wave function -^exp{^ j^adi]. After such a transform the magnetic potential 
A is replaced by A + c J^{da/dq) dt, but the electromagnetic tensor Fjk (1.2b) remains 
unchanged. So, without loss of generality one can assume that 

a = 0, E = -c~ 1 dA/dt, B = curl A 

The quantum dynamical equations have the following form: 
- in the nonrelativistic case 

dip 1 
ih— = H(q,p(t))ip, where H ~ —g jk (q)p jPk ; (8.1a) 
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- in the relativistic case 

h 2 ^ + H 2 (q,p(t))iP = 0, where H 2 ~ c 2 (g jk (q)p jPk + m 2 c 2 ). (8.1b) 

In the latter case the metric gi is assumed to be non-negative definite; the symbols ~ in 
(8.1 a,b) mean that some terms of order h, h 2 could be added to the Hamilton function 



52, M 20 



In equations (8.1) for each fixed time t the operators p(t) = —ihd/dq — -A(t,q) and 
q = q satisfy relations (1.3) with time-dependent tensor 

F jk (t, q) = e -e ki iB\t, q), q G M 3 , j, k = 1, 2, 3. 

The time derivative of the operators p(t) in (8.1) is the following: 

-p j {t) = eE j {t A ), j = 1,2,3. (8.2) 

So we see that the electric field is responsible for "dynamical evolution" of the quantum 
magnetic algebra (1.3). 

Let us introduce two-point electric potential 

H 

P(t,q,q')= / E(t,q)dq (8.3) 

Jq 

(the integral is taken along the straight line segment), and also the two-point magnetic 
potential 



1 r 

«(^g,?') = l JT \q-q\B(t,q)dq. 

\Q — 1 I Jq' 



U) 



We stress that these potentials are different from those used by Valatin [I7| in the time 
dependent case, since in our present definitions there is no integration over the time 
variable. Time and space are separated because we study the Cauchy problem in time. 

Lemma 8 The relation holds: 

83 Ida „ , 

Now from (8.2) and from composition formulas (2.11) we obtain the following state- 
ment. 

Proposition 6 (i) The time derivative of any Weyl function in quantum coordinates q, 
pit) is given by 

-ihj t f(q,p{t)) =f e (q,p(t),t), 

where 

f e (q,p,t) =eP(t,L q ,Rq)f(q,p), 
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and L q = q + \ihd p , R q = q — \ihd p are operators of the regular representation (2.10). 
(ii) The composition of two Weyl functions is given by 

[f2(q,p(t))} ■ [fi(q,p(t))\ = HQ,P(t)), k = f 2 (L q ,L p (t))fi 
where L p (t) = p — \ihd q — ^a(t, L q , R q ). 

The solution of the evolution problem (8.1a) has the general form 

^(t,q)=U t (q,p(t))U 



t=o 



In view of Proposition 6 equations for symbol U l are the following 

- ih^ + e(3(t, L q , R q ) + H(L q , L p {t)]\u\x) = 0, U° = 1. 



(8.5) 



(8.6) 



The operator acting on U l can be represented (as in (7.5)) via a symbol Tin over R t x 
x M.y. In the same way as in (7.6) we have 

H h = H + 0(h 2 ), Ho(t,x,y) = e(3(t,l q (x,y),r q (x,y)) + H(l q (x,y),l p (t,x,y)), 

(8.7) 

where x=(q,p),y= (y q ,y p ), and 

J<? = r q = q + ±y p , l p = p + \y q - \a{t,l q ,r q ). 

(Of course, here we just re-state identities (5.4) in a new notation.) As in Section 7, the 
WKB-solution of (8.6) has the form (7.7), (7.11), where (X,Y) is now the trajectory of 
the Hamiltonian system 



X = 



dH 

dy 



(t,X,Y), Y = 



dx 



(t,X,Y), X 



= x 



Y 



t=o 



= 0. 



t=o 



(8.8) 



The solution of problem (8.1b) with additional Cauchy data dip/dt\ = can also 
be constructed in the form (8.5), where U l satisfies the equations 



f) 2 

in— - e /3(t, L q , R q )\ U f = H\L q , L p {t))U\ 



d_ 
dt 



t=0 



The WKB-approximation has the form 

1 

± 



Z^ = ^exp{i5 ± }« t ± + 0(ft), 



(8.9) 



where the phases S± and amplitudes u l ± correspond (by formulas (7.11)) to the Hamilton 
function H, of type (8.7) with ± signs in the definition of H. The Hamiltonian system 
(8.8) again plays the basic role. 

The difference in Hamiltonian system (8.8) from the earlier (7.10) is that function Ho 
in (8.8) now depends on l q , l p and on r q as well. So, r p is not an integral of motion for (8.8). 
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Thus instead of dynamical system (7.13) we get now two systems: one for 7 = l(t, X, Y) 
and another for A = r(t, X, Y). They are the following: 

OH 

% = ^(l q ,l P ), (8.10) 

OH e 
7p = —q^^v lp) ~ ~B(t, 7 9 ) x % + eE(t, 7,), 



and 



A 9 = 0, (8.11) 
A p = e£(t, A g ) 



t=o 



= x°. 



t=0 



with one and the same initial condition 7 

The function H(q,p) has the following form: H(q,p) = g^ k (q)pjPk/2m in the nonrela- 
tivistic case and H(q,p) = ±Ca/ g^ k (q)PjPk + rn 2 c 2 in the relativistic case. 

Note that (8.10) is the standard dynamical system for charged massive particle in the 
electromagnetic field. The additional system (8.11) can be interpreted as the dynamical 
system for a particle of charge e and mass m = 00. The appearance of this additional 
"virtual particle" is due to the presence of the electric field E. 

The phase S of the WKB-solution is given by (7.11); hence, 



S(t,x)= [ YXdt-e [ /3(t, 7 J,A*)df- / H(^)dt. (8.12) 
Jo Jo Jo 



Lemma 9 The following identity holds: 

f YXdt= f uj f = [ oj + - [ B(t,q)dqAdq, (8.13) 

where the membrane E* C IR 6 = T*IR 3 is constructed from the two trajectories 7 = 7* 
(8.10) and A = A* (8.11) and from the magnetic wing with vertices [7*,x, A*], where x = 
(q,p) (see Fig. 15, right picture). The projection of E* onto IR 3 £/ie membrane E* 
constructed by a piece of the trajectory {7q|0 < t < t} and by the chord [?°,7 g ] u^^/i 
middle point q. 

This is the membrane area interpretation of the first term in (8.12). The second term 
in view of definition (8.3) can also be written as two-dimensional area, but in extended 
space-time: 

-el f3dt = e dt I E(t,q)dq = e EdqNdt. (8.14) 

Jo Jo Jq° J ~ 



E 



t,q° 
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Here S ti9 o is a membrane in M f x whose boundary consists of the trajectory (the world 
line) {(t, 7*) | < t < t} + the chord [7*, q°] + the straight time-segment {(t, q°) | < i < 
t}. 




Figure 17. 

Now one can combine (8.13), (8.14) and apply the Stokes theorem to transform the 
integration area to be of the most elegant geometry. Let us denote by S tiX the membrane in 
Mt x ]R® whose boundary is constructed by the world line of the given particle {(t, 7*) | < 
t < t}, the world line of the virtual infinitely heavy particle {{t, A*) | < t < t}, and also 
by the magnetic wing with vertices [A*, x, 7*] (see Fig. 17). We refer to £ 4>x as a dynamical 
membrane. 

Proposition 7 The WKB-phase of symbol U l in (8.5) can be represented as 

S{t,x)= [ (coo + F) - [ #( 7 f ) dt, (8.15) 

Js^ JO 

where ujq = ^Jdx A dx, the 2-form F is given by (1.2b), £ t)!E is the dynamical membrane 
in M7 = M. t x WL®, and 7* is the solution of classical dynamical system (8.10). 



Remark 5 The closed 2-form iop = Uq + F, which appeared in (8.15), generates a contact 



structure on 



l t XJ 



54, 55]. The "virtual" system (8.11) is the characteristic system 



for up. More precisely, the vector field on R corresponding to (8.11) is 
This is the null- field for up'- 



v \uj f = 0, 
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and the flow of vo preserves ujf'- 

C VQ LJ F = 0. 

Here we denote by C the Lie derivative and use the sign J for the contraction of a vector 
field and a form: v\uj (u) = u(u,v) for all u. If one denotes by v H the vector field on M 7 
corresponding to (8.10) 

then 

vh\ojf = dH — Vq(H) dt, C Vh uj f = d(v (H)) A dt. 

Here Vo(H) = eEdH/dp; so we see how the electric field E determines the "nonconserva- 
tion" properties of the charged particle dynamics in the contact space R 7 = R t x R^- 

Now let us return to the WKB-representation (7.7), (7.11) of symbol U l and calculate 
the Jacobian J = T>X/T>x° . The trajectory X of system (8.9) is given now by a modifi- 
cation of (7.12): X(x°,t) = |(7*(a; ) + A*(x )) +v t (x°), where v* is the same as in (7.12). 
Since the solution A = A* of (8.11) is easily calculated: A^ = q°, A* = p° + e J* E(t, q°) dt, 
we derive 



J = det 

Here 



2 v- i - 1 ; i i Qt D t 



(8.16) 



c <9g° 



The function a 5 is determined by a s {t,q,q') = \{ct{t,q,q') + a(t,q',q)), where the two- 
point magnetic potential a is given by (8.4). The point x° everywhere in these formulas 
has to be expressed via t, x by means of the equation 

x = i(7*(x ) + AV)) + v\x°). (8.17) 

This equation is uniquely solvable while the Jacobian is positive 

J>S>0, te[0,T]. (8.18) 

So, we conclude with the following result. 

Theorem 8 The symbol U l of the evolution operator (8.5) solving the equation of motion 
(8.1a) or (8.1b) can be represented {for sufficiently small time (8.18)) in the WKB-form 
(7.7) or (8.9) over the contact space R 7 = R t x R^. The phases S are given by membrane 
formula (8.15) and amplitudes u l = J~ x l 2 by (8.16). 
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Remark 6 Of course, the contact space R can be symplectified (see |5l, E5l) up to 
R 8 = (R po © E t ) x with symplectic form u' F = d Po A dt + l F . The dynamkal membrane 
H t ,x in (8.15) can be blown up to a membrane E£ x in such a way that the path 7* is put 
on the level p = —H, and the path A* is put on the level p = 0; so, the summand f* H dt 
in (8.15) is included into the membrane area, and altogether one obtains 



S(t,x) = / Up. 

JT,' 
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